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We examine supergravity models with grand unification at Mq possessing 
R parity invariance. Current data has begun to significantly constrain the 
parameter space. Thus for mSUGRA, accelerator data places a lower bound 
on mi/2 of mi/2 ~ 300 GeV while astronomical data on the amount of relic 
dark matter narrowly determines mo in terms of mi/2 (for fixed value of tan (3 
and Aq) due to co-annihilation effects. Additional new data could fix the 
parameters further. Thus the parameter space is sensitive to the muon mag- 
netic moment anomaly, 5a^, and if (5a^ lies la above its current central value, 
it would exclude mSUGRA, while if it lies la below (but is still positive) it 
pushes the SUSY spectrum into the TeV domain. The Bg decay 
is seen to be accessible to the Tevatron Runll with branching ratio sensitiv- 
ity of Br[Bs l^^fJ' ] > 6.5 x 10^^ with 15fb~^/detector, and a value of 
7(14) X 10^® obtainable with 2fb^^ would be sufficient to exclude mSUGRA 
for tan/? < 50(55). Measurements of Bg f^'^tJ'~ can cover the full mSUGRA 
parameter space for tan/5 > 40 if 6a^ > 11 x 10~^°, and combined mea- 
surements of Bg and (or alternately the gluino mass) would 
effectively determine the mSUGRA parameters for /i > 0. Detector cross sec- 
tions are then within the range of planned future dark matter experiments. 
Non-universal models are also discussed, and it is seen that detector cross 
sections there can be much larger, and can be in the DAMA data region. 
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1 Introduction 



It is generally expected that the Standard Model will break down at energies above LEP, 
and signals of new physics will occur for energies ~ 100 GeV - 1 TeV. The nature of this 
new physics is one of the crucial question of particle physics. Simultaneously, astronomical 
data has determined with good accuracy the amount of dark matter in the universe, 
though the nature of that dark matter remains one of the crucial questions of astronomy. 
Supersymmetric theories with R parity invariance offer an explanation to both puzzles as 
well as a window on the cosmology of the very early universe at times t pa 10~^ sec. 

Unfortunately, supersymmetric models depend upon a large number of parameters, 
and even the simplest model, mSUGRA, depends on four parameters and one sign. But 
fortunately, supersymmetry applies to a wide number of phenomena, and it is now becom- 
ing possible to significantly restrict the parameter space. The general MSSM, with over 
100 free parameters (63 real parameters) is not very predictive. We consider here there- 
fore, models, based on supergravity grand unification at Mq ~ 2 x 10^^ GeV (which have 
both theoretical and experimental motivation). We examine first the current status of the 
simplest model, mSUGRA |]l], and what might be obtained from future measurements 
of the muon magnetic moment, — 2, the Higgs mass, and the B decay Bg — > 
We will then look at non-universal models with non-universality in the Higgs or third 
generation of squarks and sleptons. For all these cases, the lightest neutralino, Xi, is the 
dark matter candidate, and this will also strongly constrain the SUGRA parameter space. 

2 mSUGRA Model 

We briefly review the mSUGRA model which depends on four parameters and one sign, 
and thus is the most predictive of the SUSY models. We take these parameters to be the 
following: mo (the universal scalar soft breaking mass at M^); mi/2 (the universal gaugino 
mass at Mg); Aq (the universal cubic soft breaking mass at Mg); and tan (3 =< H2 > 
/ < Hi > (the ratio of the two Higgs VEVs at the electroweak scale). The sign of the 
Higgs mixing parameter fi (appearing in the superpotential as ^HiH2 ) is the remaining 
parameter. (We note at the electroweak scale that the Xi and lightest chargino, xf , masses 
are related to mi/2 by m^o = 0.4mi/2 and m-± = 0.8mi/2-) We examine this model with 
the following parameter ranges: tjiq < 1 TeV; mi/2 < 1 TeV; 2 < tan/3 < 55; and 
Aq < 4mi/2- The above bound on mi/2 corresponds to a gluino mass range of mg < 2.5 
TeV, which is also the upper mass reach for the LHC. 

In the early universe, the neutralino can annihilate via s-channel Z, and H, and A 
neutral Higgs bosons [h is the light Higgs, and H{A) are the heavy CP even (odd) Higgs), 
and also through t-channel sfermion diagrams. However, if a second particle becomes 
nearly degenerate with the Xi? oi^s must include it in the early universe annihilation 
processes. This leads to the co-annihilation phenomena. In SUGRA models with gaugino 
grand unification, this accidental near degeneracy occurs naturally for the light stau, fi. 



1 



One can see this analytically for low and intermediate tan/3, where the renormalization 
group equations (RGE) can be solved analytically. One finds for cr, the right selectron 
and the Xi at the electroweak scale the results: 

mg^ = m^ + {)Ahm\i2 — siTn?0wM'^cos2(3 (1) 
m^o = 0.16m^/2 

where the last term of Eq(l) is approximately (37GeV)^. Thus for ttiq = 0, becomes 
degenrerate with Xi at mi/2 =370 GeV, and co- annihilation thus begins at mi/2 =(350GeV 
- 400GeV). As mi/2 increases, mo must be raised in lock step (to keep 
More precisely, it is the fi which is the lightest slepton and this particle dominates the 
co-annihilation phenomena. In general, co-annihilation implies that one ends up with 
relatively narrow allowed corridors in the mo — mi/2 plane with mo closely correlated with 
mi/2, increasing as mi/2 increases. 

Dark matter detection of Milky Way neutralinos incident on the Earth depends upon 
the neutralino - proton cross section. For detectors with nuclear targets containing heavy 
nuclei, this is dominated by the spin independent cross section. The basic quark diagrams 
involve s-channel squarks, and t-channel h and H diagrams. Thus o"-j^o_p decreases with 
increasing mi/2 and mo (which as we have seen above increase together), and also increases 
with tan [3 (due to the Higgs couplings to the d quark). Thus the maximum cross section 
will occur at high tan/?, and low mi/2, mo. 

In order to carry out calculations in SUGRA models accurately, it is necessary to take 
into account a number of corrections, and we list the important ones here: 

We use two loop gauge and one loop Yukawa RGE from Mq to the electroweak scale 
(which we take as (^1^2)^''^); and QCD RGE below for the light quark contributions. Two 
loop and pole mass corrections are included in calculation of m/^. One loop corrections to 
irii, and m^- Q are included, which are important for large tan (3. Large tan (3 NLO SUSY 
corrections to 6 ^ 57 [§, ^ are included. All stau-neutralino co-annihilation channels are 
included in the relic density calculation with analysis valid for the large tan/? regime 

Note that we do not include Yukawa unification or proton decay constraints in the 
analysis as these depend sensitively on post-GUT physics, about which little is known. 
In fact in string or M-theory analyses with grand unification, while the unification of the 
coupling constants occur as in SUGRA models, the Yukawa unification or proton decay 
constraints do not generally hold [p!0|] . 

3 Current Experimental Constraints 

In order to see what the currently allowed parameter space is, one must impose all the 
present experimental constraints. However, three of these acting together produce signif- 
icant limitations, and we mention these here: 
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Figure 1: mo — mi/2 plane showing the allowed parameter region for tan/5 = 10, Aq =0, 
fi > 0. The short vertical lines give cr^o_p = 5 x 10^^ pb (left line) and 1 x 10^^ pb (right 
line) . 



(1) Higgs mass. The current LEP bound on the light Higgs is rrih > 113.5 GeV |TT 
However, the theoretical calculation of m/j may have a 2 -3 GeV error, and so we will 
conservatively interpret this bound to mean (theory) > 111 GeV. 

(2) b ^ s'-f decay. There is some model dependence in extracting the b s'j branching 
ratio from the CLEO data, and so we will take a relatively broad range around the current 
CLEO central value [|T3: 



1.8 X 10"* < Br{B X,7) < 4.5 x 10"* (2) 

(3) Xi relic density. The relic density is measured in terms of = p/pc where p is 
the mass density, pc = SH"^ /SnGjy and H = {100km/ sMpc)h is the Hubble constant. 
Analyses of the cosmic microwave background now gives a fairly accurate measurement 
of the amount of CDM, i. e. ficDM/^^ = 0.139 ± 0.026 0. We take a 2a range around 
the central value: 

0.07 < fi^o/i^ < 0.21 (3) 

These three constraints now combine to greatly restrict the mSUGRA parameter space. 
Thus the rrih bound for low tan (3 and the b s'j constraint for higher tan (3 produce a 
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Figure 2: mo — rni/2 plane showing the allowed parameter region for tan/? = 50, Aq = 0, 
fi > O.The short vertical lines give (from left to right) cr^o_p = 5 x 10^® pb, 1 x 10^^ pb, 
2 X 10-9 pb. 



lower bound on mi/2 over the entire parameter space of 

mi/2 ~ (300 - 400) GeV (4) 

and consequently m^o ~ (120 — 160) GeV. This means that most of the parameter space 
is in the fi — Xi co-annihilation domain in the relic density calculation, and thus to satisfy 
the relic density bound , tjiq is approximately determined by mi/2 (for fixed tan/5, Aq). 
This implies that as mi/2 increases, so does mo, and so generally on has that (Jj^o_p is a 
deceasing function of mi/2. 

We consider first /i > 0. Figs. 1 and 2 exhibit the effects discussed above. Thus in 
Fig. 1 for tan/3 = 10, Aq = 0, one sees that the Higgs mass bound requires mi/2 300 
GeV, and one sees the narrow mo band allowed by the co-annihilation effects. The short 
vertical lines show the expected dark matter detection cross sections of cr^i-p = 5"^ pb 
(left) and 1 x 10"^ pb (right). Thus mo is determined by mi/2 to within ~ 40GeV. Fig. 
2 shows the corresponding situation for tan/3 =50. Here the 6 — > 57 produces the low 
energy cut off of mi/2 ~ 350 GeV. There is also a bulge in the dark matter allowed domain 
at low mi/2 where co-annihilation has not yet set in. The cross sections are (left to right) 
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5x10 ^pb, 1x10 ®pb and 2 x 10 ^ pb. These cross sections are within range of future 
planned detector experiments. 
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Figure 3: Minimum values of function of mi/2 for /x < 0, tan/? = 6 (dashed), 

8 (dotted), 10 (sohd) and 25 (large dash) [7]. 



For /i < an accidental cancellation in cr^o_p occurs over a wide range of tan/5 |^ 
giving large regions where cT-^o_p < 10 
detectors. This is exhibited in Fig. 3 [ 
as one scans the allowed parameter space, for tan/5 =6 (dashed), 8 (dotted), 10 (solid), 
25 (large dash). In this case the spin independent cross section can fall below the very 



pb and hence probably inaccessible to future 
I, where the minimum cross sections are plotted 



small spin dependent cross section where such cancellations do not occur |16 



There exist now two other experiments that might restrict the parameter space even 



more: The BNL E821 ||T^ fi'/^ — 2 experiment measuring the muon magnetic moment, and 
the decay Bs which may be observable at the Tevatron RunllB (or possibly in 

B-factories), and we turn to consider these next. 



4 Muon magnetic Moment Anomaly 

The BNL E821 experiment has measured the muon magnetic moment with exceedingly 
high accuracy. When compared with the calculations of expected from the Standard 
Model (with corrected sign of the hadronic scattering of light by light contribution 
there remains a small discrepancy: 
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6a^ = 27(16) X 10 



-10 



(5) 



which is a 1.6 a effect. While this is not enough to presume that a real effect has been 
discovered, it is still interesting to examine the effects such an anomaly would have for 
two reasons: first the errors will shortly be significantly reduced, and second SUGRA 
models imply the existence of such an anomaly of just this size. 
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Figure 4: Allowed parameter space in the mo -mi/2 plane for tan/3 = 40, Aq = 0, 
yU > 0. The upper (blue) region would be eliminated if < 11 x 10^^°, and the line for 
= 27 X 10"^° is shown. 



Much of the uncertainty in the calculation of comes from the hadronic part, a^(had). 
However, it is possible to express this by a dispersion relation in terms of integrals over the 
experimental cross section cXe+e- had. Recent experiments from from CMD-2, VEPP- 



2M and Beijing ||T9|, ^ ^ have re-measured these cross sections with greatly improved 
accuracy. Further, the BNL E821 experiment has about six times more data which they 
are currently analyzing which should greatly reduce the statistical part of the error. Thus 
one may expect the error in Eq.(5) to be reduced by a factor of two or more in the near 
future. If the central value of the anomaly were to remain unchanged, the effect would 
become more statistically significant. 

From the theoretical side, it has been known for some time that mSUGRA predicts 



an important contribution to [^2|, This is illustrated in Fig. 4 which shows the 
allowed regions in the mo - mi/2 plane and assumes that the entire anomaly is due to 
SUSY. The upper right (blue) region corresponds to less than la below the current 
central value, while the diagonal line corresponds to the central value, which falls at the 
lower edge of the allowed part of the parameter space. The current central value thus 
corresponds to a relatively low mass SUSY spectra, easily accessible to the LHC. In fact, 
too big a value of delta (i. e. delta ~ 40 x 10~^°) would be sufficient to exclude 
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Figure 5: Illustrated 95% CL. limits on the branching ratio for Bg fJ'^fJ'^ at CDF in 
Run II as a function of integrated luminosity. Solid (Case A) and dashed (Case B) curves 
are based on different assumptions on the signal selection efficiency and the background 
rejection power [30]. 



mSUGRA |2J]. On the other hand, an less than la below the current central value 



would imply a heavy SUSY mass spectrum. Further, if Sa^ is positive, then > pq, |26[|, 
and this would eliminate the very low dark matter detection cross sections shown in Fig. 3. 
Thus the final value of 6a^ will play an important role in deciphering the SUSY parameter 
space. 



5 B. 



11^ 11' 



The Bg — i> fi'^fi^ decay offers an additional window for investigating the mSUGRA pa- 
rameter space. This process has been examined within the MSSM framework ^ 
and more recently using mSUGRA E^. We consider here predictions for this decay 
for mSUGRA, but include all the current experimental constraints listed in Sec. 3 
(which are necessary to see what predictions occur|5^). The Bg ^~^fi~ decay is 



of interest since the Standard Model prediction for the branching ratio is quite small 
2^{Br[Bg /i+yu"] = (3.1 ± 1.4) x 10"^), while the SUSY contribution can become 
quite large for large tan/5. This is because the leading diagrams grow as (tan/5)^ and 
hence the branching ratio as (tan/3)^. What further makes this decay interesting is that 
it is possible to find a set of cuts so that CDF (and probably also DO) may be able to 
observe it in Run2B (with 15 fb~^ data). Fig. 5 shows the CDF limit on the branching 
fraction as a function of the luminosity BO]. (The solid curve is a conservative estimate. 
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Figure 6: Branching ratio for Bs f^^fJ'^ (three dashed hnes from left to right: 7 x 10~®, 
2 X 10^*^, 1 X 10^*^) for tan/9 = 50 in the mo-mi/2 plane. Other mSUGRA parameters are 
fixed to be y4o = and fi > 0. The three short solid lines indicate the o"-j^o_p values (from 
left: 0.05 xlO^^ pb, 0.004 xlO"^ pb, 0.002 xlO"^ pb). The vertical dotted lines label 
Higgs masses [30]. 



and the dotted curve a more optimistic possibility.) One sees that CDF will be sensi- 
tive to a branching ratio of Br > 1.2 x 10~^ (and the combined CDF and DO data to 
Br > 6.5 X 10^^). mSUGRA analysis then shows that CDF would be sensitive to this 
decay for tan (3 ~ 30. 

We can now examine the effects of the combination of all data. Figs. 6 and 7 [show 
the parameter space for tan [3 = 50 and tan (3 = 40 respectively for Aq = and /i > 0. 
One sees that there is an upper bound on the Br[Bs Ai"*"/^"] to be consistent with 
mSUGRA, and a branching ratio > 7(14) x 10^*^ would be able to exclude mSUGRA 
for tan (3 < 50(55). As can be seen from Fig. 5, such a branching ratio could be seen 
with only 2 fb~^. More generally, Fig. 6 shows that the entire parameter space can be 
covered for tan/3 = 50, = if > 11 x 10~^°, and from Fig. 7 also for tan/? = 40 
using the combined CDF and DO data. The effect of varying Aq is shown in Fig. 8 where 
Aq = — 2mi/2, and tan/5 = 40. Here again the full parameter space can be covered. 

One sees from the above graphs that future measurements should be able to determine 
the basic parameters of the mSUGRA model. Thus when the lines for 5a^, Br[Bs 
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Figure 7: Branching ratio for Bg fJ'^fJ'^ (three dashed hnes from left to right: 1.9 x 10^^, 
1 X 10^®, 0.7 X 10^^) at tan/3 = 40 in the mQ-mi/2 plane. Other mSUGRA parameters 
are fixed to he Aq = and /i > 0. The three short solid lines indicate the cr^o_p values 
(from left: 0.03 xlO'^ pb, 0.002 xlO"^ pb, 0.001 xlO'^ pb). The vertical dotted lines 
label Higgs masses [30]. 



and TTih intersect the allowed dark matter band at a single "point" for a given Aq, 
this will determine mo (approximately), mi/2 and tan/3. (It may be better to use the 
gluino mass in place of rrih since the parameter space is very sensitive to m/j.) Predictions 
of the SUSY spectrum and dark matter detection rates would then follow. 



6 Non-Universal Models 

One can generalize mSUGRA by allowing non-universal soft breaking at Mq in the third 
generation of squarks and sleptons and also in the Higgs sector. If universality of the gaug- 
ino masses is maintained, then stau-neutralino co-annihilation will still play an important 
role. However, new effects can occur, since the non-universality effects the size of the fi 
parameter. The fi parameter governs the Higgsino content of the Xi and as fi"^ decreases 
(increases), the Higgsino content increases (decreases). Since (Jj^o_p depends on the inter- 
ference between the Higgsino and gaugino parts of Xi, c"x°-p ^^^^ correspondingly increase 
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Figure 8: Branching ratio for Bg fJ^^fi at tan/5 = 40 in the mo-mi/2 plane for 
Aq = — 2mi/2 and /i > 0. The two short sohd hnes indicate the cr^o-p values (from left: 
0.005x10-6 pb, 0.001 xlQ-^ pb). The vertical dotted lines label Higgs masses[30]. 



(decrease). A second effect also occurs. As the Higgsino content of the Xi increases, then 
the Xi - Xi - ^ coupling is strengthened, allowing a new annihilation channel to become 
important (in addition to the fi-x? co-annihilation channel). As a simple example we 



consider the case where at Mq one chooses ttih^ 



rrir 



'1 + 52), and all other soft breaking 



masses universal. Fig. 9 shows the allowed region in the itlq - mi/2 plane for tan/3 = 40, 
62 = 1, Ao = 0,fi>0. One sees the usual narrow stau-neutralino co-annihilation band 
a relatively low mg, but in addition there is a higher tuq (and low mi/2) region satisfying 

all constraints due to the new Z-channel annihilation process. Fig. 10 shows (T-j^o p as a 

function of mi/2 for tan/3 = 40, Aq = mi/2, fi > 0. The Z-channel corridor now reaches 
up to the DAMA data region |]^ for low mi/2 {ttl^o ^ 0.4mi/2), and so if the DAMA 
results are confirmed, it would point to a non-universality of this type. 



7 Conclusions 

We have considered here SUGRA models with R parity invariance and grand unification 
at Mg ~ 2 X 10^6 GeV. Current data has begun to constraint these models significantly. 



10 



1000 



800 



> 

O 600 



400 



200 



\ 



400 600 800 



1000 



Figure 9: Allowed regions in the mo - r7ii/2 plane for 82 = 1, tan/5 = 40, Aq = 0, fi > 0. 
The lower band is the ususal fi-Xi co-annihilation channel, and the upper region is the 
new Z boson s-channel annihilation region. 



For mSUGRA (and many non-universal models), accelerator data {rrih and b sj) place 
lower bounds on mi/2 such that mi/2 ~ SOOGeV (or m^o ~ 120GeV), while astronomical 
data on the amount of dark matter narrowly determines mO in terms of mi/2 (for each 
Aq and tan/3). 

Thus additional new data could begin to fix the parameters of mSUGRA parameters, 
and so determine the mass spectrum expected at the LHC, dark matter detection rates, 
etc. In particular, the muon magnetic moment anomaly and the branching ratio 
Br[Bs /i/i] could combine with rrih or a measurement of rrig to fix the parameters 
completely. A more accurate determination of should be available shortly, and we have 
seen that the Tevatron Runll should be sensitive to the Bg decay for Br > 1.2(0.65) x 10^^ 
for the CDF (or combined CDF and DO) data, and this would cover a large part of the 
parameter space for tan (3 > 30. Non-universal models in the Higgs and third generation 
soft breaking masses offer the possibility of a new Z boson s-channel annihilation in the 
early universe. Such possibilities can significantly increase the neutralino-proton detection 
cross section up the DAMA region of 10^^ pb. 
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Figure 10: o"^o-p as a function of mi/2 {m^o = OAmi/2) for tan/3 = 40, ^ > 0, rrih > 114 
GeV, Aq = mi/2 for ^2 = 1. The lower curve is for the fi — Xi co-annihilation channel, 
and the dashed band is for the Z s-channel annihilation allowed by non-universal soft 
breaking. The curves terminate at low mi/2 due to the b ^ s'j constraint. The vertical 
lines show the termination at high mi/2 for 6a^ > 11 x 10~^'^[24]. 
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